Proton radioactivity within a generalized liquid drop model 
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The proton radioactivity half-lives of spherical proton emitters are investigated theoretically. 
The potential barriers preventing the emission of proton are determined in the quasimolecular 
shape path within a generalized liquid drop model (GLDM) including the proximity effects between 
nuclei in a neck and the mass and charge asymmetry. The penetrability is calculated in the WKB 
approximation. The spectroscopic factor has been taken into account in half-life calculation, which is 
obtained by employing the relativistic mean field (RMF) theory combined with the BCS method with 
the force NL3. The half-lives within the GLDM are compared with the experimental data and other 
theoretical values. GLDM works quite well for spherical proton emitters when the spectroscopic 
factors are considered, indicating the necessity of introducing the spectroscopic factor and the success 
of the GLDM for proton emission. Finally, we present two formulae for proton emission half-life 
similar to the Viola-Seaborg formulae and Royer's formulae of a-decay. 
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I. INTRODUCTION 

The opportunities provided by the radioactive beam 
facilities make the study of the exotic nuclei with ex- 
treme numbers of neutrons or protons a very interesting 
topic both from the experimental and theoretical points 
of view. And studies on these exotic nuclei lead to the 
discovery of a new form of radioactivity-proton emission. 
The proton drip line represents one of the fundamental 
limits of nuclear existence and the nucleus with a very 
large excess of protons undergo spontaneous proton emis- 
sion towards stability. Besides, the rapid proton capture 
process which plays a very important role in nuclear as- 
trophysics has its inverse in the proton radioactivity from 
the nuclear ground state or low isomeric states. There- 
fore the study of the proton emission is significative. This 
proton emission from nuclei was firstly observed in an 
isomeric state of ^^Co in 1970. With the development 
of experimental facilities and radioactive beams, proton 
emissions from ground state or low isomeric states have 
been identified between Z = 51 and Z = 83 p]] and more 
proton-emitting nuclei will be observed in experiments in 
the future. 

The proton radioactivity can be used as an useful 
tool to extract nuclear structure information such as the 
shell structure and the coupling between bound and un- 
bound nuclear states [1]. Measurement on the proton 
energy, half-life and proton branching ratio (fine struc- 
ture) helps to determine the angular momentum I car- 
ried away by the emitted proton and to characterize its 
wave function inside the nucleus [H, i, g. The proton 
emission can be dealt with in the framework of WKB 
barrier penetration model since the decay process can 
be treated in a simple quantum tunneling effect through 
a potential barrier. Several approaches have been em- 
ployed to study the half-lives of spherical proton emit- 
ters, such as the distorted- wave Born approximation 
the density-dependent M3Y(DDM3Y) effective interac- 



tion [1,J0], JLM interaction fi\, and the unified fission 
model U. However, the calculations in Ref. [1, 0, Q 
did not take account the spectroscopic factor. The spec- 
troscopic factor is very important from the viewpoint of 
the nuclear structure. In this study, we calculate the 
spectroscopic factor by employing the relativistic mean 
field theory in combination with the BCS method, and 
then determine the partial proton emission half-lives of 
spherical proton emitters within the macro-microscopic 
GLDM. The calculated results are compared with the 
experimental data and other theoretical results. 

The paper is organized as follows. In Sec. II, the 
framework of GLDM and the method for spectroscopic 
factor are presented. The calculated proton emission 
half-lives are shown and discussed in Sec. III. In Sec. 
IV, two formulae for proton emission are proposed. Fi- 
nally, we present a brief summary of the present work. 



II. METHODS 

The GLDM allows to describe the processes of fusion, 
fission, light nucleus and a emission [13, [HI, [H, [H, [3, 
m, [3 [13 ■ The macroscopic energy includes the volume, 
surface, Coulomb and proximity energies: 



E — Ev + Es + Ec + Epio 



(1) 



For one-body shapes, the volume, surface and Coulomb 
energies are given by: 



Ev = -15.494(1 - l.8I^)A MeV, 



(2) 



Es = 17.9439(1 - 2.6/2)^2/3 (5/47ri?2) MeV, (3) 



Ec = 0.6e2(zVi?o)x0.5 / {V (9) / Vq) {R{d) / Rof sin 9 d9. 

(4) 
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/ is the relative neutron excess and S is the surface area of 
the one-body deformed nucleus. V{9) is the electrostatic 
potential at the surface and Vq the surface potential of 
the sphere. When the nuclei are separated : 

Ev = -15.494 [(1 - 1.8/i^)yli + (1 - 1.8/1)^2] MeV, 

(5) 



Es ^ 17.9439 



(1 - 2.6I^)AI^^ + (1 - 2.6/2)^2 



Ec = 0.6e^ZiVi?i + 0.6e^Z^/i?2 + e^ZiZa/ 



2 r^2 



MeV, 
(6) 

(7) 



The additional centrifugal energy Ei coming from the 
angular momentum of the emitted proton has been in- 
troduced: 



Ei(r) 



(8) 



Here Ai, Zi, Ri and li are the mass number, charge num- 
ber, radii and relative neutron excesses of the two nuclei. 
The relative neutron excess of a proton is fixed at /2 = 
to ensure a negative volume energy and a positive sur- 
face energy. The dimensionless quantity of I is the an- 
gular momentum carried by the emitted proton(angular 
momentum transfer), /x is the reduced mass, r is the 
distance between the mass centers. The radii Ri of the 
daughter nucleus and proton are given by plil] : 



Ri = (1.284^/^ -0.76 + 0.8Aj fm,i = 1,2. (9) 

The radii Rq of parent nucleus can be obtained with the 
volume conservation: 



l/3x 



Rq — R'l + i?2 . 



(10) 



The surface energy results from the effects of the sur- 
face tension forces in a half space. The nuclear proxim- 
ity energy iSprox has been introduced to take into account 
the additional surface effects due to the attractive nuclear 
forces between the surfaces in a neck or a gap between 
two separated fragments: 



£^Prox('') = 27 



<S>[D{r, h)/b]27rhdh, 



(11) 



where h is the distance varying from the neck radius or 
zero to the height of the neck border. D is the distance 
between the surfaces in regard and b — 0.99 fm is the 
surface width. $ is the proximity function of Feldmeier 
and the surface parameter 7 is the geometric mean 
between the surface parameters of the two nuclei or frag- 
ments: 



7 = 0.9517^(1 - 2.6/f)(l - 2.6/|) MeV-fm"''. (12) 

The partial half-life of a spherical proton emitter is cal- 
culated using the WKB barrier penetration probability. 



The decay constant of the proton emitter is simply de- 
fined as A = i^oPSp and half- life Tp = In2/A. The assault 



frequency vq has been taken as 8 x 10 
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The barrier 



penetrability P is calculated by the following formula: 



P — exp 



^j2B{r){E{r) - i;(sphere))dr 



(13) 

where i?in and -Rout are the two turning points of the 
WKB action integral. The two following approximations 
are used here: i?in ^ Ri + R2 and B{r) = /i. i?out is 
given as: 



Rout — 



ZiZ2^ 

2Q 




ZiZ2e^ 



2Q 



l{l + i)h? 
2^iQ 



(14) 



For proton radioactivity, 
given bv [sLfigj: 



the spectroscopic factor is 



Sp — Uj , 



(15) 



where uj is the probability that the spherical orbit of 
emitted proton is empty in the daughter nucleus. Fortu- 
nately, the daughter nuclei of spherical proton emitters 
are all in ground states. Thus, it is relatively easy to 
determine this spectroscopic factor by using the RMF 
theory combined with the BCS method. The RMF au- 
tomatically includes the spin-orbit interaction. It has 
received much attention due to its great success in de- 
scribing the structure of the stable nuclei [l^l, neutron- 
rich nuclei [2l[, proton-rich nucl ei p^. superdeformed 
nuclei ^23^1 and superheavy nuclei [24l. |25|. [2a | . It is now 
a standard tool in low energy nuclear structure. The 
RMF theory is well known and it will not be discussed in 
detail here. The pairing correlation is treated by the BCS 
method. We have introduced the strength of the pairing 
forces in the following forms for neutrons and protons, 
respectively [26l |: 



- |(1 
Gp = ?(1 



N 



2A 

N - Z 



) MeV 
) MeV 



(16) 

A - ' 2A ' - ■ ^1^) 
which depend on the proton number Z and neutron num- 
ber N . A is the total mass number. The NL3 parameter 
set, which has been used with enormous success in the de- 
scription of a variety of ground-state properties of spheri- 
cal, deformed and exotic nuclei [l^,!!^], is used here. Un- 
like the situation near the neutron drip line, for proton- 
rich nuclei the Coulomb barrier confines the protons in 
the interior of the nucleus. As a consequence, the effects 
of the coupling to the continuum is weaker and therefore 
for nuclei close to proton drip line, the RMF -f BCS model 
could still be considered as a reasonable approximation 
providing sufficiently accurate results [2^ . 



III. HALF-LIVES OF SPHERICAL PROTON 
EMITTERS 
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The values of angular momentum transfer /, decay en- 
ergy (5, penetrability P, spectroscopic factor Sp, the ex- 
perimental and calculated proton emission half-life are 
given in Table I. The experimental Q value, which is a 
crucial quantity to determine the decay half-life, is used 
for the calculation and the potential barrier given by the 
GLDM has been adjusted to reproduce the experimental 
Q value. 

In Ref. [6], D. N. Basu and coworkers calculated the 
half-lives of spherical proton emitters using the density- 
dependent M3Y effective interaction with the density 
of the daughter nucleus taking from phcnomcnological 
models. Madhubrate Bhattacharya and G. Gangopad- 
hyay developed this model by obtaining this density from 
mean field calculation and also used the JLM effective in- 
teraction . The results have been presented in the last 
three columns of Table I for comparison. As can be seen. 
The DDM3Y and JLM models can provide good explana- 
tion for most cases. Their is no doubt that the DDM3Y 
and JLM models are very successful because their micro- 
scopic nature includes many nuclear features. In a num- 
ber of decays, however, the results do not match so well, 
such as i-^^Tm, i^^Lu, ^^^Ta and i56Ta(isomeric state). 
For ^^^Tl(isomeric state) and ^®^Bi, the discrepancies are 
off by an order of magnitude compared with the experi- 
mental data. It is their calculations without introducing 
the spectroscopic factor that leads to a large deviation for 
some nuclei. In order to obtain more information about 
proton-radioactivity and preform more accurate investi- 
gation theoretically, we calculated the spectroscopic fac- 
tor within the RMF-t-BCS, and the results are shown in 
the fifth column of Table I. One could notice that, spec- 
troscopic factor is S'p ~ 1 at the beginning of a proton 
shell of residual daughter nucleus, but moving away, it 
decreases, becoming quite small at the end of the shell. 
This shows the shell structure plays an important role for 
the spectroscopic factor. Hence the spectroscopic factor 
includes the shell effect to a large extent. For ^^^Tl (iso- 
meric state) and ^*^Bi, the spectroscopic factors are very 
small, leading to the longer half-lives than the DDM3Y 
and JLM calculated. 

We calculate the half-lives by employing the GLDM 
taking account the spectroscopic factors. The penetra- 
tion probabilities obtained with GLDM are shown in the 
fourth column. The penetrability P stay between 10"^"^ 
and 10~^^ which are relatively very large while the range 
is narrow, compared with 10~^^ ^ 10~^^ for a-decay 
[2^ . Hence it is easy for the proton to escape from 
the proton emitter, confirming that the proton-emitting 
nuclei are weakly bound. The calculated half-lives are 
presented in the seventh column. The half-lives by the 
GLDM with the average discrepancy less than 40%, give 
better agreement with the experimental data than that 
by DDM3Y and JLM models. For nuclei ^-^^Tm, ^^°Lu, 
i^^Ta, i5'^Ta(isomeric state) and i^'^Tl(isomeric state), 
which the DDM3Y and JLM can not explain well, the 
GLDM could give a very excellent results with the de- 
viation 4%, 26%, 5%, 10% and 33% respectively. This 



indicates that the including of the spectroscopic factors 
in calculations is necessary. The quantitative agreement 
with the experimental data are better than other theo- 
retical ones which demonstrates that the GLDM with the 
proximity effects, centrifugal potential energy, the mass 
asymmetry and spectroscopic factor could be used to in- 
vestigate the proton emission successfully when the right 
Q values are given. The GLDM overestimates the half- 
life of ^^^Bi by one order of magnitude. The first reason 
possibly is the uncertainty of the Q value or the emit- 
ted proton being not in Si/2 state, and it requires further 
investigation theoretically and measurements with high 
accuracy. The second reason is, perhaps, the shell effect 
is not included in the potential barrier, though the most 
shell effect contribution has been included in the Q value 
and spectroscopic factor. It is sure that the GLDM con- 
nected with WKB approximation will quantitatively give 
more consistent results for the proton radioactivity half- 
life when the shell effect is included in the penetration 
barrier. 

Recently, the spherical proton emitter ^^^Ta was ob- 
served. Its emitted proton energy and half-life have been 
measured again (30j . With the angular momentum trans- 
fer /= 5 [3l'| and proton energy Ep = 1.444 ± 0.015 
MeV (Q = 1.453 ± 0.015 MeV), we obtain half-life of 
Tp = 3.91}'q ms using the GLDM compared with the 
experimental data of Tp — 2.9^\'i ms. The Q value is 
compatible with the half-life, indicating the experimental 
data should be reliable. The new spherical proton emit- 
ter ^^^Re was synthesized in the reaction -'^"^Gd (^^Ni, 



p4n) ^^^Re [32| and its proton emission Q value along 
with half-life have been measured recently. The Q value 
is compatible with the half-life if and only if ^ = 5 (the 
calculated value is 23lg /iS in contrast with the exper- 
imental data of 21I4 lis), which indicates the proton is 
emitted from an 7rhii/2 state agreeing with the conclu- 
sion in Ref. 1331. 



IV. NEW FORMULAS FOR PROTON 
EMISSION HALF-LIFE 



The centrifugal potential energy Ei can reduce the tun- 
neling probability and hence increases the half-life. A 
formula can be deduced to describe this relationship be- 
tween half-life and I value for proton emission, being sim- 
ilar to the formula for a-decay in Ref. [33| : 



logiorp(0-logioTp(0) + co- 



1(1 + 1) 



v/(A-l)(Z-l)A-2/3 

(18) 

The half-life Tp is measured in second, cq is slightly 
model dependent with cq = 2.5 for the GLDM. With 
the formula (|18p . we fit new formulae that could be used 
to describe half-life for proton emission. A formula is 
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proposed in the following form: 



The other formula for proton emission is given by 



logio [Tpis)] = {aZ + b)Q 



-1/2, 



logio[rp(s)] = 



Co 



V(A-l)(Z-l)A-2/3' 



(19) 



where Z and A are charge and mass numbers of the par- 
ent nucleus respectively, and Q the proton decay energy 
in MeV. The first two terms are similar to Viola-Seaborg 
formulae [H, |35[ for a-decay and the last term is exactly 
the contribution of centrifugal barrier from formula (jlSp . 
Performing a least squares fit to the half-lives of first 25 
spherical proton emitters available in Table I, we obtain 
a set of parameters for formula (fT9l) . There values are: 
a = 0.3437, b = 4.9628, c = -31.1253 and cq = 2.5950, 
with the average deviation ct=0.153 between the experi- 
mental and formula. Additionally, we obtained a set of 
parameters for half-lives of deformed proton emitters by 
employing a least squares fit to data, which include 11 
nuclei (Z = 53 — 67) that could be found in literature [l] 
or [3l|. These parameters are: a = 0.3637, b — 4.6467, 
c = —30.9299 and cq = 2.6244. The average deviation is 
W = 0.323. The half-life increases by 3'--^4 orders of mag- 
nitude when the angular momentum transfer I is changed 
from zero to five in terms of this formula. In other words, 
the half-life of proton emission is quite sensitive to the 
angular momentum I of the emitted proton, which in turn 
helps to determine the I value when half-life and Q value 
are measured. On the other hand, so many proton emit- 
ters have been observed in experiments at present due 
to the centrifugal barriers prolonging lifetimes of these 
nuclei to a great extent. 




50 55 60 65 70 75 80 85 



FIG. 1: Deviation between the formulae (|19I20|I and experi- 
mental logarithm of half-lives for proton emission. 



Co 



1(1 + 1) 
v/(A-l)(Z-l)A-2/3 



(20) 



The first three terms are similar to Royer's formulae for 
a-decay [l^ [s^. With the same method we discussed 
above, we obtained the parameter sets: a = —23.0632, 
b = -0.4225, c = 0.4170 and cq = 2.5989 for spheri- 
cal proton emitters with average deviation a = 0.183; 
a = -23.9341, b = -0.3936, c = 0.4385 and cq = 2.6167 
for deformed proton emitters with a = 0.316. The two 
formulae (|19l20p can validate each other, and get the 
more reliable results in future study. The experimental 
proton radioactivity half-lives of most nuclei can be re- 
produced within a factor of less than 2 by the above two 
formulae. The little discrepancy suggests that these two 
formulae could be used to determine the / value when 
Q and Tp values have been measured and then extract 
some useful information about nuclear structure as well 
as to calculate proton emission half-lives. The proton 
and a emission can be described by the similar formulae 
with the different parameters, which is exactly what we 
expect. Since both the two decay modes are quantum 
tunneling effect, the studies on them should be unified. 

From these two formula, it is easy to deduce the fol- 
lowing equation: 



dQ 



1 



{aZ+b)Q-^/^, or ~-cZQ 



1 



-3/2 



(21) 



which reflects the Q value dependence of half-life. By 
comparing these parameters and Q values with that in 
Viola-Seaborg and Royer's formulae for a-decay, one 
could find that the half-life is more sensitive to Q value 
for proton emission. In addition, Cn « 2.6 for proton 
emission in contrast with Co = 1.0 [33| for a-decay indi- 
cates that the centrifugal barrier is much more important 
for proton emission than that for a-decay, due to smaller 
reduced mass n (and hence the high centrifugal barrier) 
compared with that in a-decay system. These imply that 
it is quite difficult to predict the half-life of proton emis- 
sion for unknown nucleus since the Q value can not be 
obtained with a good accuracy and since the uncertainty 
on the I value is large. 



V. SUMMARY 

The proton radioactivity of spherical proton emitters 
has been analyzed in the framework of the GLDM for 
the first time. The penetration barriers are constructed 
in the quasimolecular shape path, and the penetration 
probabilities are calculated with the WKB approxima- 
tion. The penetration probabilities are relatively very 
large while its range is narrow. Therefore it is easy for 
the proton to escape from the proton emitter, confirm- 
ing the proton emitting nuclei are weakly bound. The 
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spectroscopic factors have been taken into account in 
half-lives calculations, which arc obtained by employing 
the rclativistic mean field theory combined with the BCS 
method. The spectroscopic factor is affected greatly by 
the proton shell structure, and in turn it contains shell 
effect to a large extent. For -'^^''Tl (isomeric state) and 
^*^Bi, the spectroscopic factors are very small (0.022 and 
0.011), leading to the longer half-lives than the DDM3Y 
and JLM calculated. Although the DDM3Y and JLM 
models include the appropriate considerations in the mi- 
croscopic level, the present calculations are better agree- 
ment with the experimental data than that within the 
DDM3Y and JLM models. This, indicates the consider- 
ing of spectroscopic factor is necessary for proton emis- 
sion, especially for the nuclei with residual daughter nu- 
clei at the end of the shell. Additionally, some newly 
observed proton emitters "'^^^Ta and "'^^^Rc have been an- 
alyzed. Finally, two formulae similar to Viola-Seaborg 
formulae and Royer's formulae have been proposed for 
proton radioactivity. On the one hand, they can be em- 
ployed to calculate the proton emission half-life. On the 
other hand, they can be used to determine the I value 



when the Q value and half-life have been measured and 
then extract some useful information about nuclear struc- 
ture since the decay rate is quite sensitive to I value. The 
proton and a emission can be described by the similar 
formulae with different parameters, but the half-life of 
proton radioactivity is more sensitive to Q and I values 
compared with a-decay according to our analysis, lead- 
ing to the prediction of the half-lives for proton emission 
quite difficult. 
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TABLE I: Comparison between experimental and calculated proton radioactivity logarithmic half-lives of spherical proton 
emitters. The asterisk(*) symbols in parent nuclei denote the isomeric states. The experimental data of ^^^Ta and ^^^Re are 
taken from Ref. fsC?] and Ref. fs^], respectively and their Q values are calculated using the measured emitted proton energies. 
Other experimental data are from Ref. 



Parent 


/ 


Q(MeV) 


Penetrability 




logiorp(s) 


logio rp(s) 


logiorp(s) 


logio^pCs) 


logiorp(s) 






expt. 






expt. 


GLDM 


DDM3Y[&| 


DDM3Y[7| 


JLM[7| 




2 


0.491 


1.280 X 10"^" 


0.999 


2.049 


1.831 


1.97 


2.27 


1.69 


^Tm 


5 


1.753 


6.759 X 10"^® 


0.580 


-5.409 


-5.656 


-5.14 


-5.20 


-5.10 


1 47m 


5 


1.071 


3.931 X 10"^^ 


0.581 


0.591 


0.572 


0.98 


0.98 


1.07 


'Tm* 


2 


1.139 


2.504 X 10"^* 


0.953 


-3.444 


-3.440 


-3.39 


-3.26 


-3.27 




5 


1.283 


3.554 X 10~^° 


0.497 


-1.180 


-1.309 


-0.58 


-0.59 


-0.49 




2 


1.317 


5.734 X lO^^'' 


0.859 


-4.523 


-4.755 


-4.38 


-4.24 


-4.24 




5 


1.255 


1.839 X 10"^° 


0.490 


-0.896 


-1.017 


-0.67 


-0.65 


-0.55 




2 


1.332 


8.262 X lO"^'' 


0.858 


-4.796 


-4.913 


-4.88 


-4.72 


-4.73 


^^Ta 


5 


1.453 


5.280 X lO"^'' 


0.422 


-2.538 


-2.410 


-4.65 


-4.67 


-4.57 


1 "ifim 

^''"Ta 


2 


1.028 


4.994 X 10~^^ 


0.761 


-0.620 


-0.642 


-0.38 


-0.22 


-0.23 


1 "iRm ■Jr 

^°Ta* 


5 


1.130 


1.793 X 10"^^ 


0.493 


0.949 


0.991 


1.66 


1.66 


1.76 


^ Ta 





0.947 


1.608 X 10"^^ 


0.797 


-0.523 


-0.170 


-0.43 


-0.21 


-0.23 




5 


1.816 


1.216 X 10"^*^ 


0.308 


-4.678 


-4.636 


~ 








2 


1.284 


2.204 X 10"^® 


0.507 


-3.046 


-3.111 


-3.00 


-2.86 


-2.87 


1 fil x~. 





1.214 


2.024 X 10"^** 


0.892 


-3.432 


-3.319 


-3.46 


-3.28 


-3.29 


^•^^Re* 


5 


1.338 


1.419 X 10"^° 


0.290 


-0.488 


-0.677 


-0.60 


-0.57 


-0.49 


1 fi4x 


5 


1.844 


7.542 X lO^^'' 


0.188 


-3.959 


-4.214 


-3.92 


-3.95 


-3.86 


165J-J.* 


5 


1.733 


1.335 X 10"^^ 


0.187 


-3.469 


-3.460 


-3.51 


-3.52 


-3.44 


1661^ 


2 


1.168 


2.624 X 10~2° 


0.415 


-0.824 


-1.099 


-1.11 


-0.96 


-0.96 


166-|-j.* 


5 


1.340 


4.887 X 10"^^ 


0.188 


-0.076 


-0.025 


0.21 


0.22 


0.30 







1.086 


1.126 X 10~^° 


0.912 


-0.959 


-1.074 


-1.27 


-1.05 


-1.07 


167 jj.* 


5 


1.261 


6.559 X 10"^^ 


0.183 


0.875 


0.858 


0.69 


0.74 


0.81 


Au 





1.469 


7.608 X 10"^^ 


0.848 


-4.770 


-4.872 


-5.02 


-4.84 


-4.86 


Au* 


5 


1.718 


4.101 X 10"^® 


0.087 


-2.654 


-2.613 


-3.03 


-3.03 


-2.96 


177^1 





1.180 


1.324 X 10"^° 


0.733 


-1.174 


-1.049 


-1.36 


-1.17 


-1.20 


177 rpj* 


5 


1.986 


1.166 X 10"^'' 


0.022 


-3.347 


-3.471 


-4.49 


-4.52 


-4.46 


ISSg; 





1.624 


1.942 X 10"" 


0.011 


-4.229 


-3.392 


-5.44 


-5.33 


-5.36 



